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Abstract. Sex-biased admixture has been observed in a wide variety of admixed populations. 
Genetic variation in sex chromosomes and ratios of quantities computed from sex chromosomes and 
autosomes have often been examined in order to infer patterns of sex-biased admixture, typically 
using statistical approaches that do not mechanistically model the complexity of a sex-specific 
history of admixture. Here, expanding on a model of Verdu & Rosenberg (2011) that did not include 
sex specificity, we develop a model that mechanistically examines sex-specific admixture histories. 
Under the model, multiple source populations contribute to an admixed population, potentially 
with their male and female contributions varying over time. In an admixed population descended 
from two source groups, we derive the moments of the distribution of the autosomal admixture 
fraction from a specific source population as a function of sex-specific introgression parameters and 
time. Considering admixture processes that are constant in time, we demonstrate that surprisingly, 
although the mean autosomal admixture fraction from a specific source population does not reveal 
a sex bias in the admixture history, the variance of autosomal admixture is informative about 
sex bias. Specifically, the long-term variance decreases as the sex bias from a contributing source 
population increases. This result can be viewed as analogous to the reduction in effective population 
size for populations with an unequal number of breeding males and females. Our approach can 
contribute to methods for inference of the history of complex sex-biased admixture processes by 
enabling consideration of the effect of sex-biased admixture on autosomal DNA. 

* Corresponding author. Phone: 650-724-5122. Fax: 650-724-5114. Email: agoldb@stanford.edu 
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Introduction 

Populations often experience sex-biased demographic processes, in which males and females contributing 
to the gene pool of a population are drawn from source groups in different proportions, owing to patterns 
of inbreeding avoidance, dispersal, and mating practices (Pusey 1987; Lawson Handley & Perrin 2007). 
In humans, sex-biased demography has had a particular effect on admixed populations, populations that 
have often been founded or influenced by periods of colonization and forced migration involving an initial 
or continuing admixture process (Mesa et al. 2000; Seielstad 2000; Wilkins & Marlowe 2006; Tremblay & 
Vezina 2010; Heyer et al. 2012). 

Genetic signatures of sex-biased admixture have been empirically investigated in a variety of human popula- 
tions. In the Americas, these include African American, Latino, and Native American populations (Bolnick 
et al. 2006; Wang et al. 2008; Stefflova et al. 2009; Tishkoff et al. 2009; Bryc et al. 2010a,b; Moreno-Estrada 
et al. 2013; Verdu et al. 2014). Sex-biased admixture and migration have also been examined in populations 
throughout Asia (Oota et al. 2001; Wen et al. 2004; Chaix et al. 2007; Segurel et al. 2008; Chaubey et 
al. 2011; Pemberton et al. 2012; Pijpe et al. 2013), Austronesia (Kayser et al. 2003, 2006, 2008; Cox et 
al. 2010; Lansing et al. 2011) and Africa (Wood et al. 2005; Tishkoff et al. 2007; Berniell-Lee et al. 2008; 
Beleza et al. 2013; Petersen et al. 2013; Verdu et al. 2013). 

Sex-specific admixture and migration processes have typically been studied using comparisons of the Y 
chromosome, which is paternally inherited, and the mitochondrial genome, inherited maternally (Seielstad 
et al. 1998; Oota et al. 2001; Wood et al. 2005; Bolnick et al. 2006; Gunnarsdottir et al. 2011; Lacan 
et al. 2011). More recently, as the Y chromosome and mitochondrial genome each represent single non- 
recombining loci that provide an incomplete genomic perspective, sex-biased admixture has been examined 
by comparisons of autosomal DNA to the X chromosome (Lind et al. 2007; Wang et al. 2008; Bryc et 
al. 2010a,b; Cox et al. 2010; Beleza et al. 2013; Verdu et al. 2013). 

The Y-mitochondrial and X-autosomal frameworks are both sensible, as both involve comparisons of two 
types of loci that follow different modes of inheritance in males and females. What has not been clear, 
however, is that autosomal data, which have not typically been viewed as the most informative loci for 
studies of sex-specific processes, can carry information about sex-biased admixture, even in the absence of a 
comparison with other components of the genome. 

We demonstrate this surprising result through an extension of a mechanistic model for the admixture history 
of a hybrid population. In a diploid autosomal framework, Verdu & Rosenberg (2011) examined contributions 
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of multiple source populations that varied through time, without considering sex specificity. Here, expanding 
on the model of Verdu & Rosenberg (2011), we develop a model that mechanistically considers sex-specific 
admixture histories in which multiple source populations contribute to the admixed population, potentially 
with varying female and male contributions across generations (Fig. 1). In an admixed population descended 
from two source populations, we derive the moments of the distribution of the fraction of autosomal admixture 
from a specific source population, as a function of sex-specific admixture parameters and time. We analyze 
the behavior of the model, considering admixture processes that are constant in time, and we show that the 
moments contain information about the sex bias. 

The model 

Several studies have described mechanistic models of admixture (Chakraborty & Weiss 1988; Long 1991; 
Ewens & Spielman 1995; Guo et al. 2005; Verdu & Rosenberg 2011; Gravel 2012; Jin et al. 2013). We 
follow the notation and style of the model of Verdu & Rosenberg (2011), studying a hybrid population, H, 
which consists of immigrant individuals from M isolated source populations and hybrid individuals who have 
ancestors from two or more source populations. The source populations are labeled S a , for a from 1 to M. 
We focus on the case of M = 2. 

We define the parameters s Qi9 _i and h g ^\ as the contributions from source populations S a and H, respec- 
tively, to the gene pool of the hybrid population H at the next generation, g. That is, for a randomly chosen 
individual at generation g, the probabilities that a randomly chosen parent of the individual derives from 
S a and H are s a ,g-i and h g —i, respectively. We define the sex-specific parameter s„ 9 _i, for 5 G {/, m}, 
as the probability that the type-5 parent of a randomly chosen individual from the hybrid population at 
generation g is from source population S a . Similarly, h s g _ 1 is the probability that the type-S parent of a 
randomly chosen individual in H at generation g is from H itself. We consider a two-sex model, using / for 
female and m for male. Thus, because each individual has one parent of each type, female and male, we 
have 



s, 



J 



,m 



i)/2 



(1) 



h. 



'9-1 = 



(^'-i + ^-i)/2. 



(2) 



The contributions to the next generation of the three source populations (Si, S2, H) sum to one: 



8l,g-l + s 2,g-l + h 



'9-1 



= 1. 



(3) 
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Similarly, the female and male contributions to the next generation separately sum to one, 

*L-i + s L-i + fcj-i - Cs-i + s la-i + K-i = L ( 4 ) 

At the first generation, g = 1, the hybrid population has not previously existed; therefore, 

h 0 = h f 0 =h o n = 0 (5) 

Si,o + s 2 ,o = s{ 0 + s{ 0 = s™ 0 + s™o = !• (6) 

The first generation has two independent parameters, s[ 0 and s™ 0 - Each subsequent generation contributes 
four additional independent parameters (s[ g _ x , 1; s™ 9 _ l5 s™g-i)i and considering the first g generations, 
there are 4g — 2 independent parameters. The model is discrete in time and assumes non-overlapping 
generations. 

Our model allows us to consider complex sex-biased admixture processes by allowing uneven sex-specific 
contributions from each source population at each generation. It reduces to the model of Verdu & Rosenberg 
(2011) when the sex-specific contributions are equal within a source population, that is, if for each g, 
s( j = s" l g _ 1 and 823-1 = s 2g-i- We perform similar computations to those of Verdu & Rosenberg 
(2011), illustrating that in certain cases, our results reduce to those obtained when sex specificity is not 
considered. 

We let L be a random variable indicating the source populations of the parents of a random individual 
from the hybrid population, H. L takes its values from the set of all possible ordered parental combinations, 
{Si Si, S\H, S1S2, HSi, HH, HS2, S2S1, S2H, S2S2}, listing the female parent first. We assume random mat- 
ing in the hybrid population at each generation, so that the probability that an offspring has a particular 
pair of source populations for his or her parents is simply the product of the probabilities for having the 
female and male parents (Table 1). 

We define the fraction of admixture, the random variable H ag ^, as the probability that an autosomal genetic 
locus in a random individual of sex 5 from the hybrid population in generation g ultimately originates from 
source population a. The sex-specific fractions of admixture are related to the total fraction of admixture 
H a ,g from source population a in generation g by H a g — {H a _ g j + i? Q , 3 , m )/2. 

Under the model, we derive expressions for the moments of the fraction of admixture. Autosomal DNA 
is inherited non-sex-specifically and from both parents; therefore, female and male offspring have identical 
distributions of admixture, and H a g t and H a g m are identically distributed. Each of these quantities 
depends on both the female fraction and the male fraction of admixture in the previous generation, but 

4 



Downloaded from http://biorxiv.org/on September 18, 2014 



conditional on the previous generation (that is, on H ag -ij and £T Q , s _i )m ), they are independent. For 
our two-population model, we consider the non-sex-specific fraction of admixture, H\. g ^, treating 8 here as 
representing either / or m, (but retaining the same meaning throughout). The quantity -ffi, s ,<5 depends on 
both sex-specific fractions of admixture from the previous generation, H\^ g -\j and ifi i9 _i, m . 

Distribution of the admixture fraction from a specific source 



The definition of the model parameters and the values from Table 1 allow us to write a recursion relation 
for the fraction of admixture from source population 1 for a random individual of sex S from the hybrid 
population at generation g, or H± tg ^. For the first generation, g = 1, we have 



1 


if £ 


= Si Si, with P[L 


— Si Si] 


_ „/ a m 
— s l, 0*1,0 


1 

2 


if £ 


= Si S 2 , with P[L 


= S1S2} 


„/ „m 

— *1,0*2,0 


1 

2 


if £ 


= S 2 Sx, with P[L 


— S2S1] 


„/ „m 

— *2,0*1,0 


0 


if £ 


= S2S2, with P[L 


= S2S2] 


„/ „m 

— a 2,0*2,0 



For all subsequent generations, g > 2, we have 



1 


if £ 


= Si Si, with P[L 


= SrSi] 


„/ -m 

- s l,g-l s l,g-l 


l+fli.B-l.m 

2 


if £ 


= SiH, with P[L 


= S x iJ] 


- s l,g-l n g-l 


1 

2 


if £ 


= Si S 2 , with P[L 


= S1S2 


„/ „m 

- s l,g-l s 2,g-l 


1 + ^1,9-1,/ 

2 


if £ 


= if Si, with P[L 


= PSi] 


- n g-l s l,g-l 


Hl,g — 1,/ + H"l,g — 
2 


if £ 


= HH, with P[£ = 


= ffff] = 


- h$ h m 

- n g-l n g~l 


2 


if £ 


= H S 2 , with P[£ 


= HS 2 ] 


— hf i m 

— n g-l b 2,g-l 


1 

2 


if £ 


= S 2 S X , with £[£ 


= S2S1 


„/ e m 

- s 2,g-l s l,g-l 


Hl,g-l,m 
2 


if £ 


= S 2 i£ with P[£ 


= S 2 P] 


— h m 

— b 2,g-l' l g-l 


0 


if £ 


= S2S2, with P[L 


= S2S2 


„/ „m 

— s 2,g-l*2,g-l 



Using eqs. (7) and (8), we can analyze the distribution of the fraction of admixture as a function of the time 
g and the parameters s{ g , s™ g , g , s™ g . Under our model, -ffi.g.a takes its values in Q g = {0, 1/2 9 , 1 — 
1/2 S , 1}. Therefore, using eqs. (7) and (8), and recalling that -Hi,g,/ and -ffi lffj m are identically distributed, 
for a value q in the set Q g , we can compute the probability P(£fi. 9 .5 = 9) that a random individual from the 
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hybrid population at generation g has admixture fraction q. For g — 1, Q\ = {0, h, 1}, and 



s{ fi sf fi if g = 1 

P(ffi,i, 4 = <z) = «j s ( oS - + s ( qS - if q = I 
0 / 0 m 



For all subsequent generations, # > 2, for g in Q 5 , 



' h m V 






_ x + h 


+ (4, 


_ x + h 



n,g-l,Si 



23- 



if q = 0. 



23q-r 
29- 1 



(9) 



(10) 



The function I s is defined for all values of q in Q g , and is equal to 



s l.g-l s l,g-l 



if Q = 1 



i 9 (q) = { 



= / „m i „/ „m :f „ _ 1 

, 1.<7-I 6 2,n-1 ' *2, o-l A l,o-l 11 1 — 2 



5 2,o-l' 5 2, 0 -l 



0 



if q = 0 
otherwise. 



(11) 



In eq. (10), we calculate the probability distribution of i?i lff ,5 by taking a sum over all possible parental 
pairings at the previous generation that would lead to an admixture fraction q at generation g. Only three 
values of q allow for a history without a single hybrid ancestor — q = 0, q = |, and q = 1 — producing 
the terms in eq. (11). When there is no sex bias and s[ g _ 1 = sT.g-ii ^l-i = ^<^-i> an d s 2. 9 -i = s 2™ g -i> 
eqs. (9)-(ll) reduce to the corresponding eqs. (3)-(5) from Verdu & Rosenberg (2011). 

Eqs. (9)-(ll) can be used to analyze the behavior of the distribution of i?i i£( _i,5 over time. In Figure 2, 
we consider constant admixture processes after the founding of the hybrid population (s s a = s s a for each 
a E {1,2}, S 6 {/, 77i }, and g > 1), plotting ¥(Hi y9y s) for the first six generations, as computed recursively 
using eq. (10). In Figure 2A and 2B, we consider a hybrid population founded with equal contributions from 
source populations Si and S2, but with no further contributions after g = 1. In both of these cases, the 
distribution of the autosomal admixture fraction contracts around the mean of | . However, whereas Figure 
2A has equal contributions from each sex in the founding generation, Figure 2B has a large initial sex bias. 
We see that the width of the distribution is smaller with the sex-biased contributions, despite equality of 
the total contributions si 0 and s 2 .o- 
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In Figures 2C-2E, we consider admixture scenarios in which the founding of the hybrid population is followed 
by constant contributions from the source populations over time, si =0.1 and S2 = 0.3. Because the two 
source populations contribute after the founding, the distribution does not contract around the mean as in 
Figures 2A and 2B. Also, because the total contributions from Si and S2 are unequal, the distribution of 
Hi g s is no longer symmetrical. Rather, because the contribution from S2 is greater, the distribution is 
shifted toward zero. 

Figures 2C and 2D have the same continuing contributions for g > 2, with no sex bias in the founding 
generation for Figure 2C, and a large initial sex bias for Figure 2D. Despite different founding contributions, 
Figures 2C and 2D have similar distributions of H\ g ^ after a few generations. In Figure 2E, the hybrid 
population is founded without a sex bias and with equal contributions from the two source populations. 
The total contributions s\ and S2 are the same as in Figures 2C and 2D, but unlike in Figures 2C and 2D, 
the continuing contributions are sex-biased, with s{ ^ s™ and s{ 7^ s ™- Even when s\ and S2 are held 
constant, the distribution of -Hi.g.a depends on the s a . Notably, the probability of -ffi,6,<5 = 0 drops from 
0.157 in Figure 2C to 0.000 in Figure 2E. Similarly, the probability of -ffi,6,<5 = 1 drops to zero in Figure 2E 
as well. With these reductions at the extremes, we see a rise in the probability of intermediate values for 



Expectation of the fraction of admixture 

Using the law of total expectation, we write the expectation of the fraction of admixture from source 
population 1 for a random individual of sex 5 in population H at generation g as a function of conditional 
expectations for all possible pairs of parents L, 





(12) 



Si Si 



SiH 



S1S2 



HSi 



HH 



HS 2 



S2S1 



S 2 H 



S2S2 
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We can simplify this recursion relation. For g = 1, 

E[fli,i,*] - F{L = 5i5i)E[ifi,i, 4 |L = Si Si] 

+ P(L = SiS 2 )E[Fi >v |L = SiS 2 ] 
+ V(L = S 2 Si)-E[H 1>1)B \L = S a Si] 
+ F(L = S 2 S 2 )E[JJi )1>5 |i = S 2 S 2 ]. 

For all subsequent generations, g > 2, we have 

E[J? 1)fli(5 ] = P(L = SiSi)E[£Ti ifl , a |i = Si Si] 

+ P(L = S lJ ff)E[Fi lfl , a |L = SiH] 
+ P(£ = SiS 2 )E[Fi >5!(5 |L = SiS 2 ] 
+ P(L = HSi)E[H ltgtS \L = HSt] 
+ ¥(L = HH)E[Hi <g , s \L = HH] 
= HS 2 )E[H 1 , g , s \L = HS 2 ] 
= S 2 S 1 )E[H Ug , s \L = S 2 Si] 
+ V(L = S 2 H)E[H ltgiS \L = S 2 H] 
+ ¥(L = S 2 S 2 )E[H hg . 5 \L = S 2 S 2 ]. 



Using eqs. (7) and (8), for the first generation, g = 1, we have 

E[Hi,m] = *£o<oE[l] 

+ ( s l,0 s 2^0 + s 2,O s ™o) 

+ 4. 0 4'; 0 e[o]. 

For all subsequent generations, g > 2, we have 



E 



„/ „m 



1 + Hi„-\. 



K^sT^E 



1 + ffi 





T 


) E 






2 



#1,0-1,/ + Hl q - 



g-l,m 



H 



l,g— l,m 



Hi. 



9-1,/ 



L,-A-iE[o]. 
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Recalling eqs. (1), (2), and (4), we can simplify the expectation of the fraction of admixture in a random 
individual of sex d from the hybrid population. For g — 1, eq. (15) gives 

E[fli,i,tf] - s h0 , (17) 

the same expression found by Verdu & Rosenberg (2011, eq. 10). For g > 2, by eq. (16), 

E[ifi, 9)5 ] = S i, g -i + - [h^ElH^j] + h^nHt^J) . (18) 

Because H\ g j and H\ s ^ m are identically distributed, recalling eq. (2), we can simplify the expectation using 
E[i?i, ff ,/] = E[-ffi, g .m] = where (5 is left as an unspecified sex (/ or to). For g > 2, the expectation 

of the fraction of admixture from source population 1 is 

E[fr w ] - si, 9 _! + /i 9 _iE[ir 1>9 _ 1)5 ]. (19) 

We see in eqs. (17) and (19) that the expectation of the fraction of admixture for a random individual of 
sex 5 from the hybrid population at generation g, E[Hi^ g ^], depends on the total contributions of the source 
populations (Si, 5*2, H) at each generation, si, g _i and h g —\, and not on the sex-specific parameters, s{ g _ 1 , 
s , ™ g _ 1 , hg_ x , and h g n _ 1 . This recursion (eqs. (17) and (19)) is the same as in the non-sex-specific model of 
Verdu & Rosenberg (2011, eqs. 10 and 11). 



Higher moments of the fraction of admixture 



We can write a general recursion for the higher moments of the fraction of admixture from population Si in 
a randomly chosen individual of sex 5 from the hybrid population. For k > 1, in the first generation, g — 1, 
we have 



irk 

-"1,1,5 



l fc if L = 


Si Si, with P[L = 


Si Si] 


~ *i,cri,o 


(§)* ifi = 


Si S 2 , with P[L = 


S1S2] 


— „/ „m 

— *1, 0*2,0 


ifL = 


S 2 Si, with P[L = 


S2S1] 


— „/ „m 

— *2, 0*1,0 


0 fe if L = 


S2S2,with P[L = 


S2S2] 


— „/ „m 

— *2, 0*2,0 



(20) 
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For all subsequent generations, g > 2, 



rrfe 

H hg,s 



(§)' 



1.3-l./+ g l,g-l, 

2 



(§)* 

/ gl,3-l, 

V 2 



if L 


= Si Si, with P[L 


— *SxiSi 


- s f 

~ 1,9 


-1 1,5- 


if L 


= Si if, with P[L 


= Si#] 


- 


l n g-l 


if L 


= Si S 2 , with P[L 


= S1IS2. 
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= HSi] 
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~ S 2,9 


m 

-1 S 1,9- 


if L 


= S 2 H, with P[L 


= S 2 H] 




1^9-1 


if i 


= S2S2, with P[L 


= S2S2 




m 

-l S 2,g- 



(21) 



As in the case of k = 1, we use the law of total expectation to write a recursion for higher moments of the 
distribution of the fraction of admixture for all k > 1. Using the values for the recursion for the fraction of 
admixture, eqs. (7) and (8), for the first generation, g = 1, we have 



J a m , / m 
'1,0*2,0 ' *2,0*1,0 



s{ 0 s^ 0 E[Q h 



(22) 



For g > 2, we have 



E [ H LA =4,9-l S ™9-l E [ lfe ] 



1 + Hi, a - 
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2~ 



' t o-l 6 2,g-l JI1 



4,g-l4':g-lE[0 fc ]. 



1 + #l,g-l,/ 



g l,9-l,/ 

2 



(23) 



Recalling eq. (3) and noting that /iq = 0, we use the binomial theorem to simplify the recursion for the 
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moments of Hi_ g j. For g = 1, we have 



E [Hl hS ] = s{ fi sT fi + 



2 k 



(24) 



For g > 2, we have 



WlH k 1 - J e m i a 1.9-l°2, ff -l -r °2,g-l°l,g-l 
fi L-"l,9,<5j — s l,g-l s l,g-l 



„/ (.in 
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2 fc 
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2 fe 
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iVi 




2 fe 
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I./ _m 



(25) 



Because H\ g j and H\^ g m have the same distribution, we can simplify the fcth moment of the distribution 
of the fraction of admixture from Si, for S £ {/, m}, to give 



F \H k 1 - J « m I S L-l S 2!g-l + 5 2 ,g-l s l,-. i 
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2 fc 







(26) 



For fc = 1, eqs. (24) and (26) should produce the expectation that we have already derived for k = 1. For 
k = 1, using eqs. (1), (2), and (4), eq. (24) gives 



/ m I f m 
Tw-r-rj i f m , 1,0 2,0 2,0 1,0 



(27) 



which matches eq. (17). For g > 2 and = 1, eq. (26) gives 

Tff\fJ 1 _ J e m | S l,g-l S 2,g-l + S 2 ,g-l S l,g-l + S l,g-l ft g-l + 

^[■"l.s.aj — *l,g-l*l,g-l "T 2 

/ in x a x -f -f «„_iSi,„_i T «n-l s 2,g-l + s 2,g-l n q-l 



E[H 1)B - llS ], (28) 
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which simplifies to match eq. (19). Finally, with equal contributions in each population from females and 
males, so that s[ g _ 1 = s™ g _ 1 = Si )S _i and 4,g-i = s ™ y g-i — s 2, g -i, eqs. (24) and (26) reduce to eqs. 16 
and 17 from Verdu & Rosenberg (2011). 

Variance of the fraction of admixture 

When k = 2, eqs. (24) and (26) produce a recursion for the second moment of H\ tg s- Recalling eqs. (l)-(6), 
for g = 1, we have 

r 9 i 4 ol 1 + s ?o) + s'l'V 1 + 4 o) , x 

E [Hf^g] = — ^ ^ (29) 



For g > 2, we have 



E \H( 



2 ] _ J c m , s l,g-l°2,g-l T ^>2,g-l a l,g-l 
3 ,<5j - s l, 9 -l s l,s-l + 4 

, ''o-l 6 l,g-l 



(l + 2E[ J ff 1)fl _ 1)/ ]+E[ J ff 1 Vx, / ]) 

E [ H l a -i,f] + 2E[ffi, s -i,/]E[^i, ff -i,m] +E [^L_i, m ]) 



4 

n g-\ n Q-\ 



4 

+ ^=f^E [if x Vx, m ] + /1 -^E [fl?^ Xi/ ] . (30) 
Recalling that Hi g j and Hi gm are identically distributed, eq. (30) simplifies to give 

E [Hi s ] = Sl <3-^ 1 + S ^9-^ + ^g-li 1 + s l, g -l) 



4 

+ ^i+ii'Vi E[fli8 _ M] + (E [i?i iff _i i 5]) 2 

J./ I Lin 

Using the definition of the variance V [H hg>s ] = ®[Hi,g, s ] - {HH^gj}) 2 , and eqs. (17), (19), (29) and (31), 
for the first generation, for the variance of the fraction of admixture, we have 

4 n(l - 4 n) + sfo(l - sfo) 
V[H ltU ] = — h ° (32) 

For all subsequent generations, g > 2, we have 

v[ff 1 S l,g-l( 1 ~ flfl-l) + S hg-l( l - S l)g-l) 4,g-l hf g-l + ^-l^-l pn, 1 

h l i (1 - hi ,) + /i™ i(l - K i) 2^1+^1 
+ -*=^ 9 — -. 9 ^ g — (E [H hg _ hS ]) 2 + 9 ' 1 A 9 — Y[H hg _ 1<5 ]. (33) 
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With no sex bias, so that s{ = s™ g = Si t9 and g — s™ 9 = S2. g , eqs. (32) and (33) are equivalent to eqs. 22 
and 23 from Verdu & Rosenberg (2011). 

The recursion for the variance of the fraction of admixture of a random individual of sex S from the hybrid 
population is dependent on the variance from the previous generation, the expectation from the previous 
generation, and its square. By contrast with the expectation, the variance of the fraction of admixture 
depends on the sex-specific contributions from the source populations. 

Eqs. (32) and (33) are invariant with respect to an exchange of all variables corresponding to males (super- 
script m) with those corresponding with females (superscript /). Thus, while the variance is affected by the 
sex-specific admixture contributions, it does not identify the direction of the bias. Despite the dependence of 
the variance of the autosomal fraction of admixture on sex-specific contributions, under the model, the sym- 
metry demonstrates that autosomal DNA alone does not identify which sex contributes more to the hybrid 
population from a given source population. This result is reasonable given the non-sex-specific inheritance 
pattern of autosomal DNA. 

Special case: a single admixture event 

Using the recursions in eqs. (17), (19), (32), and (33), we can study specific cases in which the contributions 
are specified. We first consider the case in which the source populations 5*1 and 5*2 do not contribute to 
the hybrid population after its founding: s[ g = s" l g = „ — s™ g = 0, and h — [h^ g + h g n )/2 = 1, for all 
g > 1. As before, at the first generation, the hybrid population is not yet formed, and ho = 0. Therefore, 
Si,o + s 2 ,o = s{ 0 + s(o = s '™o + 4™o = 1- 

Under this scenario, we can derive the exact expectation and variance of the autosomal fraction of admixture 
of a random individual from the hybrid population. In the case of a single admixture event, the expectation 
of the admixture fraction is equal to the expectation at the first generation, because the further contributions 
are all zero. Using eq. (19), si. 9 -i = S2.g-i = 0 for all g > 2. Therefore, from eq. (17), in the case of a single 
admixture event, for all g > 1, 

E[H 1>gtS ] = 5i,o- (34) 

The expectation of the autosomal fraction of admixture from source population Si is constant over time, 
and it depends on the total — not the sex-specific — contribution from the source population Si. As in the 
general case in eq. (19), for a single admixture event, a sex bias does not affect the expectation. Because 
the source populations provide no further contributions after the founding generation, unlike in the general 
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case, the mean admixture fraction does not change with time. 

Using eqs. (32) and (33), because s{ = s™ = — s™ = 0 for all g > 2, the variance of the fraction of 



For s{ o = s™0' by eqs. (1) and (2), the variance matches eq. 25 of Verdu & Rosenberg (2011). 

With a single admixture event, the variance decreases monotonically, and its limit is zero for all parameter 
values. Individuals from the hybrid population only mate within the population, decreasing the variance by 
a factor of two each generation. Thus, eq. (35) predicts that the distribution of the admixture fraction for 
a random individual in the hybrid population contracts around the mean, converging to a constant equal to 
the mean admixture from the first generation. 

In eq. (35), considering all possible pairs (s{ 0 , s^g), with each entry in [0, 1], the maximal W[Hi ig< s] occurs 
at (s{ 0 , s™ 0 ) = (| > f)) a scenario with equal contributions from the two source populations, and no sex bias. 
At the maximum, the variance is V[i?i. s . §] — l/2 9+2 . Four minima occur, at (s{ 0 ,s™o) = (0,0), (0,1), (1,0), 
and (1, 1), cases in which all individuals in generation g — 1 have the same pair of source populations for 
their two parents, and in later generations, all individuals continue to have the same value of Hi^ g $. In these 
cases, V[i?i iSi 5] = 0. 

Figure 3 plots the variance in eq. (35) as a function of the sex-specific parameters s\ 0 and s™ 0 for three 
values of g. For g = 1, a maximum of W[Hi,u] = 1/8 occurs at (s{oi s ™o) = (I)!); an< ^ a minimum, 
V[-^M s] — 0, at (s{ 0 , s™ 0 ) = (0, 0), (0, 1), (1, 0), or (1, 1). After one generation of mixing within the hybrid 
population, with no further contributions from the source populations, the maximum and minima occur 
at the same values of (sio: s ™o)' but the variance is halved (Fig. 3B). That is, for a given set of values 



0»i,o.*i5>). V[ff li2 , 4 ] = V[ff liM ]/2. Similarly, for g = 8 in Figure 3C, V[JT 1)8 ,a] = V[i^ M ]/2 7 . By g = 8, 



the hybrid population is quite homogeneous in admixture, and the variance of the admixture fraction has 
decreased to near zero for all sets of founding parameters. Therefore, the admixture fraction distribution is 
close to constant, with Hi 8 s ~ Si,o- 

We can analyze the dependence of the variance on the sex-specific parameters by considering constant 
total contributions si.o and allowing the sex-specific contributions to vary, constrained by eq. (1) so that 
0 < s{ 0 ,s" 1 q < min(l, 2si : o)- Rewriting eq. (35) in terms of Si ; o and s{ 0 , 



admixture follows a geometric sequence with ratio |. For all generations g > 1, 




(35) 



29+ 1 




(36) 
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From this expression, it is possible to observe that given a constant Si,o in [0,1], the maximal variance is 
produced when s{ 0 = s'™ 0 = Si t Q. The minimal variance occurs when (s{ 0 ,s™ 0 ) = (0, 251,0) or (251,0,0) for 
s i,o < jj or (l,2si.o — 1) or (2si,o — 1, 1) for si,o > \- This minimum only takes the value V[7?i i£ ,,5] = 0 
when 51,0 equals 0, g, or 1. 

For the specific case of si Q = 1, the total contribution for which the maximal variance occurs in Figure 3, 
we illustrate the variance at several locations in the allowed range for s{ 0 and s™o (Fig- 4). Four scenarios 
are plotted with the same total founding contribution from source population 1, sx,o — 5, but with different 
levels of sex bias. As the female and male contributions become increasingly different, the initial variance 
decreases. The largest variance for si ( o = \ occurs at s{ 0 = s™ 0 = si.o, with no sex bias. The minimum 
occurs when males all come from one source population and females all from the other. In this extreme 
sex-biased case, the variance is zero constantly over time, as each individual has a male parent from one 
population and a female parent from the other, and an admixture fraction of \. 

Special case: constant nonzero contributions 

Next, we consider the case in which an initial admixture event founds the hybrid population, and is then 
followed by constant nonzero contributions from the source populations. After the founding, for each g > 1, 
all admixture parameters are constant in time: s s a g — s s a for each a € {1,2} and 5 G {/, m}, and h s g = h s 
for each S. Thus, we have parameter values for the founding, and constant continuing admixture parameters 
s{, s™, s 2 j ancl S T- Each parameter takes its value in [0, 1], as do si and S2- By contrast, h takes its value 
in (0, 1). The case of h = 1 is a single admixture event, analyzed above. The h = 0 case is trivial because 
the hybrid population is re-founded at each generation, and the distribution of the admixture fraction thus 
depends only on the contribution in the previous generation. Therefore, we require S1+S2 7^ 0 and S1+S2 7^ 1- 
Individually, however, h* and h m can each vary in [0, 1], as long as they are not both zero or one. 

The recursion for the expectation of the autosomal fraction of admixture, eqs. (17) and (19), is equivalent 
to that derived by Verdu & Rosenberg (2011). Therefore, the closed form of the expectation is equivalent as 
well. From Verdu & Rosenberg (2011) eq. 30, we have 

, Si,o, .9 = 1 

m ltg>6 ] = { (37) 

siyi^ + si 1 ^-, g>2. 
We can use the same method as Verdu & Rosenberg (2011) to simplify the second moment. Under the 
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special case of constant contributions across generations, for g = 1, eq. (29) gives 

<o(l + ^ 0 ) + ^o(l + <o) 



mii, s ] 



For g > 2, eq. (31) gives 



TO- \tt2 1 . S {(l + ^)+^"(l + ^) , ^fe m + fe^r F r ff 1 
hfh m o /)/ -L. fr" 1 

+ ^- (E [H^Af + ^4^ E [^L-M] • 



(38) 



(39) 



Because this equation is a non-homogenous first-order recurrence with the form 



(40) 



we can use Theorem 3.1.2 of Cull et al. (2005) to solve for a unique solution for E[H? $], as in Verdu & 



Rosenberg (2011). For the initial condition, we have 

a 0 = E[H 11S ] = — — . 

We define A = (h f + h m ) /4 = h/2, and for all g > 2, we have 
sf (1 + sf ) + sf (1 + s{) s{/i m + /i-^sj 



hfh ni 

E[ff liff _ M ] + __ (E[ff li9 _ li4 ]) 2 . 



4 2 1 y ' J 2 

Using the expected admixture fraction from eq. (37), we can simplify eq. (42). For all g > 2, 



1 - /l 9 " 
1 - h 



Therefore, using Theorem 3.1.2 of Cull et al. (2005), we have a unique solution for E[H? s ]: 

a 0 , g = 1 



a a (2) 



E 



4 



( s l,0 



M" 2 -I- o-, ! 



(41) 



(42) 



(43) 



(44) 
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Eq. (44) can be simplified by separating the sum and summing the resulting geometric series: 



ata, 



2/ 



= 1 

> 2, 



(45) 



where cxq is defined in eq. (41), and 

S! + s{s™ S!{s{h m + h* sf) 



At 
A 2 
A 3 

A 4 



sfh f h m 



(2-h) (l-h)(2-h) (l-h) 2 (2-h) ] 



h f h r - 



Q'O 



■Si 



1 - h 
(s{h m + h f s^ 



si 



h V 1 - h 

S\h 



h f h n ' 
2h 2 



1 - h 
si / si 



81,0 



- Sink - 



Si 



(l-h)(2-h) 
si 



(l-h)(2-h) 



8i + s{a? 
2-h ' 



1 - h V 1 - h 



2sin 



(46) 
(47) 
(48) 

(49) 



When s{ = s™ and s{, — s™, A\, A%, A3, and A4 are equal to the corresponding quantities in eqs. 39-42 
of Verdu & Rosenberg (2011). Therefore, without sex bias, the closed form of the second moment of the 
admixture fraction, eq. 45, is equal to eq. 38 in Verdu & Rosenberg (2011). 

Using the relation V[H ltg>s ] = E[Hl g S ]-(E[H hgiS ]) 2 and eqs. (37) and (45), for the variance of the autosomal 
fraction of admixture, we have 



4 ' 

\Hi, g ,s] = { r g-i 

A 1 +A 2 h9- 1 + A 3 + A 4 J2( 2h Y 

i=l 



3=1 



(50) 



> 2. 



For h ^ h , we have 



s lo( i - s {,o)+ s "o( i - s ro) 



V[Hl,g, S ] 



A 1 + A 2 hs 



4 

-1 



A, 



3 +^( ? rf ! )](r-(^- i + ^) 2 , 



(51) 



<?>2. 



For h = |, eq. (50) gives 



a i,o( 1 - a i,o)+ s ro( 1 ~' s i!o) 



V[lfi, fl ,«] 



A! + A 2 (i)«-i + ^ + ^ _ 1} _ [ 2si + (si o _ 2si) ( i )S -i 

Eqs. (50)-(52) simplify to eqs. 43-45 of Verdu & Rosenberg (2011) when s{ — s™ and S2 — s-j 



.9=1 
3>2. 



(52) 
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Limiting variance of admixture over time 

Figure 5 illustrates the variance of the autosomal fraction of admixture as a function of g when the con- 
tributions from the source populations are constant over time, computed using eq. (50). The figure shows 
that if the continuing contributions are held constant, then the long-term limiting variance does not depend 
on the founding parameters. Unlike in the hybrid isolation case, under the scenario of constant, nonzero 
contributions from the source populations over time, h ^ 0 and h ^ 1, a nonzero limit is reached. Applying 
eq. (50), we have 

\ 2 

Si 



lim V[H ltBtS ] =A ± -( —4 , (53) 

which does not depend on the founding parameters. The limit matches that of Verdu & Rosenberg (2011, 
eq. 46) in the absence of sex bias (s{ = s™ = s lf ft/ = h m — h, and s 2 = s™ = Sa). 

The maxima and minima of the limiting variance 

Using eqs. (1), (2), (4) and (46), the limit in eq. (53) can be equivalently written in terms of the two female 
sex-specific contributions, s\ and s 2l and the total contributions from the two source populations, s\ and s 2 . 
Considering admixture scenarios with constant s%, s 2 , with s\ + s% G (0, 1], but allowing s{ and s 2 to range 
over the closed unit interval, the limiting variance depends on two independent parameters, s{^,s 2 G [0, 1], 
subject to the constraint in eq. (1): 

hm V[^ g s] = -( s ^-^£ + S is 2 ( Sl + S2 ) ^ 
g^oo 1 L ' 9 ' 01 (si+s 2 ) 2 (l + si + s 2 ) 

Treating s\ and s 2 as constants in [0,1], the critical points of eq. (54) are the same as those of 

f(s{,s f 2 ) = -( S { S2 -s f 2Sl ) 2 . (55) 

First we consider the maximum. Because f(s{,s 2 ) is always negative or zero, the maximal variance given 
si and s 2 occurs when f(s{,s 2 ) = 0, which occurs on the line s{s 2 = s 2 s\. Equivalently, recalling eq. (1), 
this line can be written as 



{ si 

7 = — = — ■ 56) 



Eq. (56) has many solutions for (s{, s™, s 2 , s™) given s\ and s 2 . One solution is {s{,s 2 ) — (si,s 2 ), which 
by eq. (1) is equivalent to (s{, s 2 ) — (s™, s 2 l ). Therefore, the limiting variance of the admixture fraction is 
maximized when there is no sex bias. Figure 6 plots two examples of the variance for constant si and s 2 , but 
increasingly different sex-specific contributions from the source populations. In both panels, the admixture 
history with no sex bias produces the greatest limit. 
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For fixed s\ and S2, however, the case without sex bias is not the only maximum of the limiting variance. 
Figure 7 plots the variance over time for four different admixture histories, each with the same total contri- 
butions si and S2, but quite different sex-specific contributions (s{, s™, s^j 8 T)- Each of the four scenarios 
plotted reaches the same limit because each provides a solution to eq. (56). Because f(s{,s 2 ) = 0, eq. (54) 
depends only on the total contributions s\ and 82- For constant S\ and S2, any admixture history whose 
contributions solve /(s{, s 2 ) = 0 has limiting variance 

lim W[Hi g s ] = -. ^ r. (57) 

L ' 9 1 ( Sl + s 2 )(i + Sl+S2 ) v > 

This maximal limiting variance depends on the total contributions from the source populations, but not on 
the sex-specific contributions; it is equivalent to eq. 47 of Verdu & Rosenberg (2011). 

Thus far, we have considered the maximal limiting variance as a function of the sex-specific parameters 
given constant total contributions s\ and Sa- We can also identify the values of s\ and S2 that maximize 
the limiting variance, considering all si,S2 G [0,1]- For each choice of Si and S2, the maximal variance 
over values of s{ and s 2 is given by eq. (57). We can therefore find the si and S2 that maximize eq. (57). 
As demonstrated by Verdu & Rosenberg (2011), given si + S2, the maximal limiting variance occurs when 
si = S2- Over the range of possible choices for si + s 2 G (0, 1), the maximum occurs when si = s 2 = \- 
Unlike in Verdu & Rosenberg (2011), however, this maximum requires the sex-specific contributions to solve 

/( s {,4) = o. 

Interestingly, one of the minima of the limiting variance occurs when S\ = s<z = |, but with /(s{, s 2 ) 7^ 0. 
Specifically, when si = s 2 = |, but all males come from one source population and all females from the other, 
(s{, s™, s 2 , s™) = (1j 0j 0j 1) or (0) 1) 1j 0)j the limiting variance in eq. (54) is zero. In this case, L — SiS 2 or 
L = S 2 Si for every individual in the hybrid population. By eq. (8), the hybrid population is founded anew 
at each generation, with each individual having admixture fraction Hi g s — \- Therefore, the population 
has zero variance. 

More generally, given s\ and s 2 , the minimal limiting variance occurs when (s{, s 2 ) = (2si, 0) or (s{, s 2 ) — 
(0, 2s2)- Given s± and S2, the limiting variance is minimized with respect to s\ and s 2 when f(s\,s 2 ) is 
smallest (eq. (54)). Because / is the negative of the square of a difference of products, it is greatest when 
one term is zero and the other is at its maximum, as at (s{, s 2 ) = (2si, 0) or (s{, s 2 ) = (0, 2S2). These points 
represent the maximal sex bias for fixed si, S2- 

If we allow s\ and S2 to vary, because a variance is bounded below by zero, any set of parameters that 
produces zero variance is a minimum. In eq. (54) if either s± = 0 or s 2 — 0, then the limiting variance of 
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the admixture fraction is zero. When only one population contributes after the founding, in the limit, all 
ancestry in the hybrid population traces to that population. 

Properties of the limiting variance 

The limiting variance of the fraction of admixture over time in eq. (53) is a function of the sex-specific 
contributions from the hybrid population, hf and h m , and source population 1, s[ and s™. Recalling 
eq. (4), the limiting variance is equivalently written as a function of the sex-specific contributions from 
source population 2, s| and s™, and either source population 1 (eq. (54)), or the hybrid population. It 
can be viewed as a function of all six sex-specific parameters (s{, s™, Sj, s™, ft/, ft m ), four of which can be 
selected while assigning the other two by the constraint from eq. (4). 

We can therefore analyze the behavior of the limiting variance as a function of two of the sex-specific 
parameters by specifying two other parameters, and allowing the final two parameters, one female and 
one male, to vary according to eq. (4). Of the four parameters we consider, using the constraint from 
eq. (4) separately in males and females, two must be male and two must be female. Because the variance is 
invariant with respect to exchanging the source populations or the sexes, the six-dimensional parameter space 
generates a number of symmetries. Figures 8-12 examine the five possible, non-redundant ways of choosing 
two populations and the corresponding male and female parameters from those populations, and holding 
two corresponding parameters fixed (either from the same sex in the two populations, or for males and 
females from one population) while allowing the other two to vary. Figure 13 then highlights an informative 
case that considers the limiting variance as a function of a male and a female parameter from different 
populations. 

Each figure shows multiple contour plots of the limiting variance as a function of two sex-specific parameters, 
for fixed values of two other parameters. Three cases plot the limiting variance as a function of the female 
and male parameters from a given population, with the female and male contributions of another population 
specified. In two other cases, parameters for a single sex from two populations are plotted, specifying the 
contributions from the other sex for those populations. 

By considering these parameter combinations, we can examine the dependence of the variance on sex-specific 
parameters and parameter interactions, as well as potential bounds on both the parameters and the variance. 
We highlight a number of symmetries in the limiting variance. The plots also illustrate the maxima and 
minima found in the previous section. 
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Properties of the limiting variance in terms of s{ and s™ 

In each panel in Figure 8, we consider the variance of the fraction of admixture as a function of s\ , the female 
contribution from Si, on the x-axis, and s™, the male contribution from Si, on the y-axis, computed using 
eq. (53). We plot the variance for fixed ft/ , the female contribution from H, and h m , the male contribution 
from H. The domain for s[ and s™ is constrained by eq. (4), with s{ taking values in [0, 1 — ft/], and s™ 
taking values in [0, 1 — ft m ]. 

The upper left plot in Figure 8 shows the variance as a function of s{ and s™, with ft/ = h m = ft = 0. 
In this setting, the hybrid population is founded anew by the source populations each generation, and s\ 
and s™ both take values from the full domain [0, 1]. For ft/ = ft m = ft = 0, the maximal limiting variance 
is lim^oo W[Hi tgj g] = |, occurring when s{ = s" 1 = s\ = 0.5. At this maximum, given eq. (4) and 
fai = — o, we have s 2 — s™ = s 2 = 0.5. Therefore, as in eq. (54), the maximal limiting variance occurs 
when female and male contributions from the source populations are equal, and the total contributions from 
the source populations are equal. 

The minima of lim^oo V[-ffi )ff)< $] = 0 occur at the four corners of the plot. At the origin, when s{ — s™ = 
si = 0, the limiting variance is zero because only S2 contributes to the hybrid population. Individuals in the 
hybrid population all have parents L = S2S2, anc ^ an admixture fraction of zero (eq. (8)). By exchanging Si 
for S2, the case of s{ = s™ = s% = 1 is similar. 

Additional minima occur at (s{, s™) = (1, 0) or (0, 1). Here, all males come from one source population and 
all females from the other. Therefore, all individuals at the next generation of the hybrid population have 
parents L = S1S2 or L — S2S1, and admixture fraction | (eq. (8)). 

For ft/ = ft™ — ft = 0, the limiting variance is symmetrical over the line s™ = s i > as a result of the symmetry 
between males and female in the variance (eq. (32)). Because the hybrid population provides no contribution 
and the variance of the fraction of admixture is symmetric with respect to source population, the variance 
is also symmetric over the lines s{ — 0.5 and s™ = 0.5. 

The columns of Figure 8 consider increasing, fixed values for , and the rows consider increasing, fixed 
values for ft m , both from {0,0.25,0.5,0.75,0.95}. All panels maintain the general shape of the limiting 
variance as a function of s{ and s™ seen for ft/ = ft™ = 0. However, as the domain for s{ and s™ shrinks 
with increasing ft/ and ft™, the location of the maximal variance changes across panels. In all cases, the 
maximum of the limiting variance occurs when s{ and s™ each lie at the midpoints of their respective 
domains, s{ = (1 - ft / )/2 and sf = (1 - ft m )/2. The magnitude of the limiting variance at each maximum 
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decreases as its location moves away from s{ = s™ = 0.5. 

For each panel, the minimum lim^oo ~V[Hi yg s] = 0 occurs when s{ and s™ are either both zero or they lie 
at the maxima of their respective domains. In these cases, only one source population contributes to the 
hybrid population, and therefore, all individuals in the hybrid population have an admixture fraction from 
Si of either 0, when s{ = s™ = 0, or 1, when s{ = 1 — ft/ and s™ = 1 — h m ((8)). The limiting variance 
is no longer zero at the two corners of each contour plot where only one of {s{, s™} is at the maximum of 
its domain; these corners, however, are minima of the variance conditional on the values of si and S2- In 
these cases, males all come from one source population and females from the other, producing a minimum 
of eq. (54) conditional on fixed si and S2 . 

As in the case of ft/ = ft™ = 0, each plot is symmetrical in reflecting over both the midpoint of the x-axis, 
s{ = (1 — ft/)/2, and that of the y-axis, s" 1 = (1 — H m )/2. The limiting variance is symmetrical with respect 
to source population (eq. (54)), and this pair of reflections corresponds to an exchange of source populations. 
For ft/ — ft'™ = 0, the line s{ — s™ generates circular contours, but as the contributions from the hybrid 
population increase, the contours become elliptical. 

In Figure 8, plots on the diagonal have equal contributions from males and females in the hybrid population, 
ft/ = h m = ft. For ft/ ^ ft m , plots above the diagonal are equivalent to those below the diagonal with an 
exchange of female for male contributions, for both s± and ft. For example, the plot with ft/ = 0.25 and 
ft™ = 0.5 is equivalent to the plot with ft/ = 0.5 and h m = 0.25 if the axes are also switched so that s™ 
appears along the x-axis and s{ is on the y-axis. 

Figure 9 plots the limiting variance as a function of s{ on the x-axis, and s™ on the y-axis, as in Figure 8, 
but we now fix values of s 2 by column and by row using eqs. (54) and (1). The maxima and minima 
occur at the same parameter values found in Figure 8, but they appear in different locations on the plots. 
For example, in Figure 9, the global maximum across panels occurs in the plot with s| = s sT = 0-5 specified, 
and s[ = s™ = 0.5. By eq. (4), this location implies — h m = 0, the plot that contains the maximal 
variance in Figure 8. In the upper left plot in Figure 9, the limiting variance is a constant zero for all s\ and 
s™, because S2 = 0 (eq. (54)). 

Whereas all panels in Figure 8 are symmetric in reflecting over the midpoints of both domains, in Figure 9, 

f f 
only the plots with s 2 = s™ are symmetric over the line s™ = s\. However, the symmetry corresponding to 

transposing males and females is visible in that a plot above the diagonal and its corresponding plot below 

the diagonal are equivalent if the axes for s{ and s™ are switched. 
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Properties of the limiting variance in terms of ft/ and h m 

Similarly to Figure 8, Figure 10 considers the limit of the variance of the fraction of admixture over time as 
a function of the four variables s{, s™, ft/, and ft™ using eq. (53). In Figure 10, each plot shows ft/ on the 
x-axis and ft™ on the y-axis, for specified values of s{ and s m , with the domains of ft/ and ft™ constrained 
by eq. (4). As in Figure 8, the plots along the diagonal are the cases with s{ = s m , and there is a symmetry 
over this line of plots in that if the values of s[ and s™ are switched, then the plots will be equivalent with 
a transposition of the axes. 

As in Figure 9, in the upper left plot in Figure 10, s{ = s™ = si — 0, and the limiting variance is a 
constant zero. In Figure 10, the maximal variance occurs at the origin (ft/ = ft™ = ft = 0) of the plot with 
s{ = s™ = s\ = 0.5. As in Figure 8, at the maximum by eq. (4), = s m — S2 = 0.5. In this case, females 
and males contribute equally. Both source populations contribute maximally to pull the distribution of the 
fraction of admixture toward the extremes of zero and one. 

Because the limiting variance is symmetrical with respect to source population, and recalling eq. (4), each 
plot in Figure 10 is equivalent to a corresponding plot in Figure 9 reflected along both the x-axis and y-axis. 
For example, the plot in Figure 10 with s{ — 0.5 and s™ = 0.25, is equivalent to the Figure 9 plot with 
si, — 0.5 and s™ = 0.25 if reflected on both the x- and y-axes. 

Figures 8-10 illustrate that the global maximum of the limiting variance occurs when the two source popu- 
lations contribute equally, the contributions from the two sexes are equal, and the hybrid population does 
not contribute to the next generation. As the parameters move from the location of the maximal limiting 
variance to the minimum, the variance monotonically decreases. 

Properties of the limiting variance in terms of s{, , and s| 

Next we plot on the x- and y-axes two parameters of the same sex from different populations. Because the 
variance is invariant with respect to transposition of females and males, we consider only females without 
loss of generality. Specifically, in Figure 11 we plot the limiting variance as a function of s[ on the x-axis 
and ft/ on the y-axis, for fixed values of s™ and ft m . In Figure 12, we plot the limiting variance as a function 
of s{ and s^, for fixed s™ and s m . For both Figure 11 and Figure 12, the domains of s{,S2, and ft/ are 
constrained by eq. (4). 

For Figure 11, the maximal limiting variance occurs in the plot with s™ = 0.5 and ft™ = 0, at (s{, ft/) = 
(0.5, 0). By eq. (4), this location is the same parameter set for the maximum in Figures 8-10. The maximum 
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in each plot occurs when (s{,/?/) — (0.5,0), but the magnitude of the variance decreases with increased 
distance from the plot with fixed s" 1 = 0.5 and h m — 0. Similarly, within each plot, the limiting variance 
decreases with distance from (s{, h?) = (0.5,0). 

In the first column of Figure 11, where s™ = 0, the line s{ = 0 produces zero variance because the hybrid 
population is homogenous, with only one source population contributing. Similarly, on the diagonal s" 1 + 
h m = 1, as st — s™ = 0 by eq. (4), the line s[ + hi =1 has minimal variance. 

In Figure 12, because of the symmetry in sex in eq. (54), the plots above and below those where s™ = s™ 
are equivalent with a transposition of axes. As in Figures 8-11, the maximal variance occurs in the plot with 
s™ = s™ = 0.5 at s{ = = 0.5. Also, similar to Figure 11, in the first column, when s™ = 0, the line 
s{ = 0 is a minimum because no contributions trace to source population Si; in the first row, when s™ = 0, 
the line s 2 = 0 is a minimum. 

Analogous to the similarity between Figures 9 and 10, by eq. (4), each plot in Figure 12 is a transformation 
of a plot in Figure 11. For example, for the plot with s™ = 0.25 and h m — 0.5 specified in Figure 11, because 
the male contributions sum to one, this panel also specifies s™ = 0.25. Therefore, we can compare this plot 
to the plot with S™ — s™ = 0.25 in Figure 12. Both show s{ on the x-axis, and using eq. (4), we can rewrite 
the y-axis in Figure 12 as = 1 — hf . 

Properties of the limiting variance in terms of non-corresponding parameters 

Finally we consider a case in which males from one population in (Si, S2, H) are compared to females from a 
different population. While multiple parameter configurations are possible, we plot one that is particularly 
informative, providing a perspective on eq. (54) beyond the observations visible in Figures 8-12. Figure 13 
plots the limit of the variance of the admixture fraction as a function of s{ on the x-axis and s™ on the 
y-axis, for fixed values of s| and s™. We rewrite eq. (54) as a function of s{,s™,S2, and using eq. (1). 
We have, 

limVLff o i = - 2 (4s? ~ 4s?) 2 + (4 + sT)(s{ + s?)(s{ + sf + s{ + sj) 

The limit depends on products of sex-specific parameters, including s{s™, as can be seen in the shape of the 
contours in Figure 13, but not in the analogous plots in Figure 9. 
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Discussion 

Our model demonstrates the potential utility of autosomal DNA in the study of sex-biased admixture 
histories. Under a framework in which admixture occurs over time, potentially with different male and 
female contributions from the source populations, we have derived recursive expressions for the expectation, 
variance, and higher moments of the fraction of autosomal admixture. For the special case of constant 
admixture over time, we have analyzed the behavior of the variance of the admixture fraction. Although 
the expectation of the autosomal admixture fraction is dependent only on the total contributions from the 
source populations, we have found that the variance of the autosomal admixture can be informative about 
sex-specific contributions. Specifically, for constant admixture over time, we have shown that the variance 
of the autosomal admixture fraction decreases as the male and female contributions become increasingly 
unequal. 

That autosomal DNA can carry a signature of sex-biased admixture might at first appear counterintuitive, 
as unlike the sex chromosomes, autosomes are carried equally in both sexes. The phenomenon can, however, 
be understood by analogy with the well-known result that increasing sex bias decreases the effective size of 
populations (Wright 1931; Crow & Dennison 1988; Caballero 1994; Hartl & Clark 2007). In a computation 
of effective size using the coalescent, for example (Nordborg & Krone 2002; Ramachandran et al. 2008), 
the sex bias causes pairs of genetic lineages to be likely to find common ancestors more recently than in a 
non-sex-biased population, as the reduced chance of a coalescence in the sex that represents a larger fraction 
of the breeding population is outweighed by the greater chance of a coalescence in the less populous sex. In 
a similar manner, if admixture is sex-biased, because lineages are more likely to travel along paths through 
populations with the larger sex-specific contributions, then the variability of genealogical paths — and hence, 
the variance of the admixture fraction — is reduced compared to the non-sex-biased case. 

Autosomal DNA, with its multitude of independent loci, potentially provides more information about the 
complex histories of hybrid populations, and the full autosomal genome might be less susceptible to selective 
pressures at individual loci than the sex chromosomes. To take advantage of autosomal information, many 
recent efforts to study sex-biased demography have focused on comparing autosomal DNA with the X chro- 
mosome (Ramachandran et al. 2004, 2008; Wilkins & Marlowe 2006; Hammer et al. 2008, 2010; Bustamante 
& Ramachandran 2009; Keinan et al. 2009; Casto et al. 2010; Emery et al. 2010; Keinan & Reich 2010; 
Labuda et al. 2010; Lambert et al. 2010; Gottipati et al. 2011; Heyer et al. 2012; Arbiza et al. 2014). Our 
study enhances the set of frameworks available for considering effects of admixture and sex bias on autosomal 
variation. 
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For a single admixture event, the expectation of the autosomal admixture fraction is constant in time and 
not dependent on sex-specific contributions. Unlike in the case of hybrid isolation, if constant nonzero 
contributions from the source populations occur over time, then the variance of the fraction of autosomal 
admixture reaches a nonzero limit, dependent on these continuing sex-specific admixture rates, but not on 
the founding contributions. In both scenarios, the variance can be informative about the magnitude of 
a sex bias in the admixture history of a hybrid population. For an arbitrary constant total contribution 
from a source population, the maximal variance occurs when there is no sex bias. The maximal variance 
across all allowable parameter values of the constant admixture model is seen when there is no sex bias, 
and equal contributions from both source populations, that is, s{ = s™ = s{ — s ™ = 0-5. Two admixture 
histories minimize the variance of the autosomal admixture fraction. First, the variance is zero when only one 
source population contributes to the hybrid population, and either all hybrid individuals have an admixture 
fraction of 0 or they all have a fraction of 1. Second, the variance is zero if all males come from one source 
population and all females come from the other source population. In this scenario, all individuals in the 
hybrid population have an admixture fraction of |. 

While the variance of the autosomal admixture fraction suggests that autosomal DNA is informative about 
sex-biased admixture, the relationship between the variance and the sex-specific parameters is complex. We 
uncovered an interesting case in which quite different sex-specific histories can lead to the same variance 
over time (Fig. 7). The variance is in fact dependent on the product of multiple sex-specific parameters, 
but not on each parameter separately (Fig. 13). In particular, when s{s™ = s^s™, we demonstrate that 
the variance is maximized (eqs. (54)- (56)). Therefore, when the equality s{s™ = s^s™ holds, the limiting 
variance depends only on the total contributions from the source populations, si and S2 (eq. (57)). The 
symmetry arises from the non sex-specific inheritance of autosomal DNA. 

We have considered two scenarios, isolation of a hybrid population after its founding, and constant contribu- 
tions from source populations to the hybrid population over time. While the admixture history of real hybrid 
populations is likely much more complex than these, our models can provide a starting point for statistical 
frameworks to estimate the parameters of mechanistic admixture models. It is noteworthy that although 
sex bias does influence autosomal variation, because autosomal DNA is not inherited sex-specifically, the 
sex that contributes more from a given source population cannot be identified with autosomal DNA alone. 
Because the X-chromosome follows a sex-specific mode of inheritance, consideration of the X-chromosome 
alongside autosomal data under the mechanistic model may be able to help differentiate between scenarios 
that produce the same variance with different choices of the sex with a greater contribution. 
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Figures 

Figure 1: Schematic of the mechanistic model of admixture over time. Two source populations, Si and 
S2, contribute both males and females to the next generation of the hybrid population H, potentially with 
time- varying proportions. The fractional contributions of the source populations and the hybrid population 
to the next generation G are Si, ff ,S2 )5 and h g , respectively. Sex-specific contributions from the populations 
are s{ , s 2 „, and s ™gj f° r females and males respectively. H a g $ represents the fraction of admixture 
from source population a G {1,2} in generation g for a random individual of sex S G {f,Tn} in population 
H. 

Figure 2: Probability distribution of the fraction of admixture from source population Si, P(Hi^gj), for 
a random individual from the hybrid population for the first six generations (eqs. (9)-(ll)). Each column 
corresponds to a specified admixture scenario, with constant contributions from the source populations over 
time after founding {s s a g = for each a G {1, 2}, 5 G {/, m}, and g > 2). 

Figure 3: The variance of the fraction of admixture, V[Hi j9i s], as a function of female and male contributions 
from source population Si in the first generation, s[ 0 and s" l 0 , in the case of hybrid isolation. (A) g = 1. (B) 
g = 2. (C) g = 8. At each generation, the variance decreases toward zero by a factor of two. Considering all 
(s{ 0 , s™o) m [0, 1] x [0, 1], the maximal variance occurs when (s{ 0 , s" l Q ) = (|, |), and the minimal variance 
occurs when (s{ 0 , s™ 0 ) = (0, 0), (0, 1), (1, 0), or (1, 1). The variance is calculated using eq. (35). 

Figure 4: The variance of the fraction of admixture, V[Hi t9t s], when contributions from the source pop- 
ulations occur only in the founding generation and the total contribution from source population 1 is held 
constant at s^o = \- The limit of the variance of the fraction of admixture over time is zero for any choice 
of (s{ n , s™o)- The magnitude of the variance, calculated from eq. (35), is inversely related to the level of sex 
bias. For all four scenarios, si = S2 = 0 and si.o = S2.0 = \- 

Figure 5: The variance of the fraction of admixture over time for constant, nonzero contributions from 
the source populations, with different levels of sex bias in the founding of the hybrid population, and 
constant, equal, and nonzero subsequent contributions from the source populations and sex for g > 1. In 
all cases, s{ = s" 1 = s{ = s™ = 0.2. The variance, calculated with eq. (50), reaches a nonzero limit 
lim^ V[Hi, gi s] = 1/24. 

Figure 6: The variance of the fraction of admixture over time for constant, nonzero contributions from the 
source populations, with different levels of sex bias, but the same total contribution from the two source 
populations. (A) s{ 0 = s™ 0 = 0.9, and s x = s 2 = 0.1, (B) s{ 0 = s\ n 0 = 0.75 and Si = s 2 = 0.2. The 
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variance reaches a nonzero limit when s{, s™, s 2 , and s™ are nonzero and constant over time. The variance 
is calculated using eq. (50). 

Figure 7: The variance of the fraction of admixture over time for constant, nonzero contributions from 
the source populations, but multiple different ratios of female to male contributions. When the sex-specific 
parameters satisfy the equation s{s™ = s™ s 2' multiple different demographic scenarios have the same 
limiting variance of the admixture fraction. In all cases, lim^oo V[i?i iS ,,5] = yg. The variance is calculated 
using eq. (33). For all scenarios s{ Q = s™ 0 = 0.95. 

Figure 8: Contour plots of the limit of the variance of the fraction of admixture over time as a function of 
s{ on the x-axis and s™ on the y-axis for specified values of ft/ by column and ft™ 1 by row. The domains of 
s{ and s™ are [0, 1 — ft/] and [0, 1 — h m ], respectively. 

Figure 9: Contour plots of the limit of the variance of the fraction of admixture over time as a function of 
s{ on the x-axis and s™ on the y-axis for specified values of s 2 by column and s" 1 by row. The domains of 
s{ and are [0, 1 — s 2 ] and [0, 1 — s" 1 ], respectively. 

Figure 10: Contour plots of the limit of the variance of the fraction of admixture over time as a function 
of ft/ on the x-axis and h m on the y-axis for specified values of s{ by column and s™ by row. The domains 
of ft/ and h m are [0, 1 — s[] and [0, 1 — s™], respectively. 

Figure 11: Contour plots of the limit of the variance of the fraction of admixture over time as a function 
of s{ on the x-axis and ft/ on the y-axis for specified values of s™ by column and h m by row. The domains 
of s{ and ft/ are bounded by the function s{ + ft/ =1. 

Figure 12: Contour plots of the limit of the variance of the fraction of admixture over time as a function 
of s{ on the x-axis and s 2 on the y-axis for specified values of s™ by column and s™ by row. The domains 
of s{ and s 2 are bounded by the function s{ + s 2 = 1. 

Figure 13: Contour plots of the limit of the variance of the fraction of admixture over time as a function 
of s{ on the x-axis and s™ on the y-axis for specified values of s 2 by column and s™ by row. The domains 
of s[ and are [0, 1 — s 2 ] and [0, 1 — s™], respectively. 

Table 1: The probabilities that an individual from the hybrid population at generation g has one of nine 
possible sets of parents from Si, S2 or H, assuming random mating. The parameter s s a g is the probability 
that the parent of sex S for a randomly chosen individual from the hybrid population, at generation g, is 
from the source population a. Similarly, the probability that this parent is from H is h s g . 
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